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Abstract. We consider discrete dynamical systems and lattice mod- 
els in statistical mechanics from the point of view of their symmetry 
groups. We describe a C program for symmetry analysis of discrete 
systems. Among other features, the program constructs and investi- 
gates phase portraits of discrete dynamical systems modulo groups of 
their symmetries, searches dynamical systems possessing specific prop- 
erties, e.g., reversibility, computes microcanonical partition functions and 
searches phase transitions in mesoscopic systems. Some computational 
results and observations are presented. In particular, we explain forma- 
tion of moving soliton-like structures similar to "spaceships" in cellular 
automata. 



1 Introduction 

Symmetry analysis of continuous systems described by ordinary or partial dif- 
ferential equations is well developed and fruitful discipline. But there is a sense 
of incompleteness of the approach since the transformations used in the sym- 
metry analysis of continuous systems — point and contact Lie, Backlund and 
Lie-Backlund, some sporadic instances of so-called non-local transformations — 
constitute negligible small part of all thinkable transformations. In this context 
finite discrete systems look more attractive since we can study all possible their 
symmetries. 

Furthermore, there are many hints from quantum mechanics and quantum 
gravity that discreteness is more suitable for describing physics at small distances 
than continuity which arises only as a logical limit in considering large collections 
of discrete structures. 

Both differential equations and cellular atomata are based on the idea of 
locality — behavior of a system as a whole is determined by interections of its 
closely situated parts. Recently 1 112 j we showed that any collection of discrete 
points taking values in finite sets possesses some kind of locality. More specif- 
ically, let us consider collection of N "points", symbolically S = {x\, ... ,Xn}- 
We call 8 domain. Each £j takes value in its own set of values Qi = { sj , . . . , sf 1 } 
or using the standard notation Qi = {0, ...,«& — 1}. Adopting Q s as symbolical 
notation for the Cartesian product Qi x • • • x Qn, we define relation on 5 as 
an arbitrary subset R s C Q s . Then we define consequence of relation R s as an 



arbitrary superset S 2 R and proper consequence as a consequence which can 
be represented in the form P a x Q s \ a , where P a is nontrivial (i.e., P a ^ (J Q ) 
relation on the proper subset a C S. We show that any relation R s allows a de- 
composition in terms of its proper consequences. This decomposition naturally 
imposes a structure of abstract simplicial complex — one of the mathematical 
abstractions of locality. Thus we call collections of discrete finite-valued points 
discrete relations on abstract simplicial complexes. 

We demonstrated also that such relations in special cases correspond to sys- 
tems of polynomial equations (if all points take values in the same set Q 
and its cardinality is a power of a prime \Q\ = p k ) and to cellular automata (if 
domain 5 allows decomposition into congruent simplices with the same relation 
on the simplices and this relation is functional) . The notion of discrete relations 
covers also discrete dynamical systems more general than cellular automata. The 
lattice models in statistical mechanics can also be included in this framework by 
considering ensembles of discrete relations on abstract simplicial complexes. 

In this paper we study dependence of behavior of discrete dynamical sys- 
tems on graphs — one-dimensional simplicial complexes — on symmetries of 
the graphs. We describe our C program for discrete symmetry analysis and re- 
sults of its application to cellular automata and mesoscopic lattice models. 

2 Symmetries of Lattices and Functions on Lattices 

Lattices. A space of discrete dynamical system will be called a lattice. Tradi- 
tionally, the word 'lattice' is often applied to some regular system of separated 
points of a continuous metric space. In many problems of applied mathematics 
and mathematical physics both metrical relations between discrete points and 
existence of underlying continuous manifold do not matter. The notion of 'adja- 
cency' for pairs of points is essential only. All problems considered in the paper 
are of this kind. Thus we define a lattice as indirected fc-regular graph r without 
loops and multiple edges whose automorphism group Aut (J 1 ) acts transitively 
on the set of vertices V(r). Sometimes we shall depict our lattices as embedded 
in some continuous spaces like spheres or tori (in this case we can talk about 'di- 
mension' of lattice). But such representations are not significant in our context 
and used only for vizualization. 

The lattices we are concerned in this paper are shown in Fig. [1] Note that 
the lattices marked in Fig. [1] as "Graphene 6x4" , "Triangular 4x6" and "Square 
5 x 5" can be closed by identifications of opposite sides of rectangles in several 
different ways. Most natural identifications form regular graphs embeddable in 
the torus and in the Klein bottle. Computation shows that the Klein bottle 
arrangement (as well as others except for embeddable in the torus) leads to 
nonhomogeneous lattices. For example, the hexagonal lattice "Graphene 6x4" 
embeddable in the Klein bottle has 16-element symmetry group and this group 
splits the set of vertices into two orbits of sizes 8 and 16. Since non-transitivity 
of points contradicts to our usual notion of space (and our definition of lattice), 
we shall not consider further such lattices. 
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Graphene 6x4 Triangular 4x6 Square 5x5 Buckyball 




Fig. 1. Examples of lattices 



It is interesting to note that the graph of hexahedron can be interpreted - 
as is clear from Fig. [2]- either as 4-gonal lattice in sphere or as 6-gonal lattice 
in torus. 



a aba 




b aba 



Fig. 2. The same graph forms 4-gonal (6 tetragons) lattice in sphere $ 2 and 
6-gonal (4 hexagons) lattice in torus T 2 . 



Computing Automorphisms. The automorphism group of graph with n ver- 
tices may have up to n\ elements. However, McKay's algorithm [4j, based on 
efficiently arranged search tree, determines the graph automorphisms by con- 
structing small number of the group generators. This number is bounded by 
n — 1, but usually it is much less. 

In Sect. [4] we discuss the connection of formation of soliton-like structures in 
discrete systems with symmetries of lattices. There we consider concrete example 
of system on square lattice. So let us describe symmetries of N x N square 
lattices in more detail. We assume that the lattice has valency 4 ("von Neumann 
neighborhood" ) or 8 ( "Moore neighborhood" ) . We assume also that the lattice 
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is closed into discrete torus Zn x Zat, if JV < oo. Otherwise the lattice is discrete 
plane Z x Z. In both von Neumann and Moore cases the symmetry group, which 
we denote by Gn x n, is the same. The group has the structure of semidirect 
product of the subgroup of translations T 2 = Zjv X Zjv (we assume Z^ = Z) 
and dihedral group D4 

G NxN = T 2 * D 4 , if AT = 3, 5, 6,..., 00. (1) 

The dihedral group D4 is, in its turn, the semidirect product D4 = Z4 x Z2. 
Here Z 4 is generated by 90° rotations, and Z 2 are reflections. The size of GnxN 
is 

|GtvxaH = 8-ZV 2 , if N^4. 
In the case N = 4 the size of the group becomes three times larger than expected 

|G 4 x 4 | = 3 x 8 x 4 2 = 384. 

This anomaly results from additional Z3 symmetry in the group G4 X 4- Now the 
translation subgroup T 2 = Z4 x Z4 is not normal and the structure of G*4 X 4 
differs essentially from ([T]) . The algorithm implemented in the computer algebra 
system GAP [S] gives the following structure 

normal closure of T 2 
G 4X 4 = ((((Z 2 x D 4 ) x Z3) x Z 3 ) x Z 2 ) xZ 2 . (2) 



Functions on Lattices. To study the symmetry properties of a system on a 
lattice r we should consider action of the group Aut (r) on the space £ = Q r of 
Q-valued functions on r, where Q = {0, . . . , q — 1} is the set of values of lattice 
vertices. We shall call the elements of £ states or (later in Sect. [5]) microstates. 

The group Aut (r) acts non-transitively on the space £ splitting this space 
into the disjoint orbits of different sizes 

N orbttB 

£= (J Oi . 

i=l 

The action of Aut (r) on £ is defined by 
where x € V(r), if (x) £ £, g G Aut (P). 

Burnside's lemma counts the total number of orbits in the state space £ 

NorMts = 1 Aut m 1 51 9^"" • 

Here N^ vcles is the number of cycles in the group element g. 
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Large symmetry group allows to represent dynamics on the lattice in more com- 
pact form. For example, the automorphism group of (graph of) icosahedron, 
dodecahedron and buckyball is SjQj and the information about behavior of any 
dynamical system on these lattices can be compressed nearly in proportion to 
|S B | = 120. 

Illustrative Data. In Table Q] we collect some quantitative information about 
the lattices from FigfT] and their automorphism groups, namely, number of ver- 
tices \V(r)\, size of automorphism group |Aut(.T)|, total number of states fi = 
\S\ = gl^( r )l (assuming q — 2) and number of group orbits N or uts m the space 
of states. 



Table 1. Lattices, groups, orbits: quantitative characteristics. 



Lattice 


m-oi 


|Aut(T)| 


Q = q\ v ™ 


■^V orbits 


Tetrahedron 


4 


24 


16 


5 


Hexahedron 


8 


48 


256 


22 


Icosahedron 


12 


120 


4096 


82 


Dodecahedron 


20 


120 


1048576 


9436 


Graphene 6x4 
Torus 


24 


48 


16777216 


355353 


Graphene 6x4 
Klein bottle 


24 


16 


16777216 


1054756 


Triangular 4x6 


24 


96 


16777216 


180070 


Square 5x5 


25 


200 


33554432 


172112 


Buckyball 


60 


120 


1152921504606846976 
^10 18 


9607679885269312 
^10 16 



3 Computer Program and Its Functionality 

We have written a C program to study different properties of deterministic and 
statistical lattice systems exploiting their symmetries. Input of the program con- 
sists of the following elements: 

— Graph F = {iVi, . . . , N n }. Ni is neighborhood of ith vertex, i.e., the set of k 
vertices adjacent to ith vertex. 

— Cellular automata branch: 

Set of local rules R — {r*i, . . . , r m }. Ti is integer number representing bits of 
ith rule. The set R includes the rules we are interested in. In particular, this 
set may contain only one rule (for detailed study). 

1 Traditionally, the icosahedral group Ih = A5 is adopted as a symmetry group for 
these polyhedra. A5 is 60-element discrete subgroup of SO (3). Adding reflections to 
A5 we get twice larger (and hence more efficient for our purposes) group S5. 
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— Statistical models branch: 
Hamiltonian of the model. 

— Some control parameters. 

The program computes the automorphism group Aut (_T) and 

— in the case of cellular automata the program constructs phase portraits of 
automata modulo Aut (_T) for all rules from R. 

Manipulating the above mentioned control parameters we can 

• select automata with specified properties, for example, reversibility, con- 
servation of a given function on dynamical trajectories, etc.; 

• search automata whose phase portraits contain specific structures, for 
example, the limit cycles of a given length, "gardens of Eden" [6] or, 
more generally, isolated cycles, " spaceships" , etc. 

— in the case of statistical lattice model the program computes the partition 
function and other characteristics of the system, searches phase transitions. 

Example of timing. 

The full run of all 136 symmetric 3-valent binary cellular automata on the dodec- 
ahedron (number of vertices = 20, order of automorphism group =120, number 
of states = 1048576, number of orbits = 9436) takes about 40 sec on a 1133MHz 
Pentium III personal computer. 

4 Deterministic Dynamical Systems 

In this section we point out a general principle of evolution of any causal dynam- 
ical system implied by its symmetry, explain formation of soliton-like structures, 
and consider some results of computing with symmetric 3-valent cellular au- 
tomata. 

Universal Property of Deterministic Evolution Induced by Symme- 
try. The splitting of the space S of functions on a lattice into the group orbits 
of different sizes imposes universal restrictions on behavior of a deterministic 
dynamical system for any law that governs evolution of the system. Namely, 
dynamical trajectories can obviously go only in the direction of non- decreasing 
sizes of orbits. In particular, periodic trajectories must lie within the orbits of the 
same size. Conceptually this restriction is an analog of the second law of ther- 
modynamics — any isolated system may only lose information in its evolution. 

Formation of Soliton-like Structures. After some lapse of time the dynam- 
ics of finite discrete system is governed by its symmetry group, that leads to 
appearance of soliton-like structures. Let us clarify the matter. Obviously phase 
portraits of the systems under consideration consist of attractors being limit 
cycles and/or isolated cycles (including limit and isolated fixed points regarded 
as cycles of period one). Now let us consider the behavior of the system which 
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has come to a cycle, no matter whether the cycle is limit or isolated. The sys- 
tem runs periodically over some sequence of equal size orbits. The same orbit 
may occur in the cycle repeatedly. For example, the isolated cycle of period 6 
in Fig. [5] — where a typical phase portrait modulo automorphisms is presented 
— passes through the sequence of orbits numbered^ as 0, 2, 4, 0, 2, 4, i.e., each 
orbit appears twice in the cycle. 

Suppose a state ip(x) of the system running over a cycle belongs to ith orbit 
at some moment to- f{x) £ 0{. At some other moment t the system appears 
again in the same orbit with the state (ft(x) = A to t (^(x)) £ Oj- Clearly, the 
evolution operator A to t can be replaced by the action of some group element 
g tot £ Aut (r) 

<p t (x) = A tot (<p(x)) = ip (gt lx) • (3) 
The element gt Q t is determined uniquely modulo subgroup 

Aut (.T; <p(x)) C Aut(r) 

fixing the state tp(x). Equation ([3]) means that the initial cofiguration (shape) 
<p(x) is completely reproduced after some movement in the space r. Such soliton- 
like structures are typical for cellular automata. They are called "spaceships" in 
the cellular automata community. 

Let us illustrate the group nature of such moving self-reproducing structures 
by the example of "glider" — one of the simplest spaceships of Conway's au- 
tomaton "Life" . This configuration moves along the diagonal of square lattice 
reproducing itself with one step diagonal shift after four steps in time. If one 
considers only translations as a symmetry group of the lattice, then, as it is 
clear from Fig. [3J (p$ is the first configuration lying in the same orbilo with tpi, 
i.e., for the translation group T 2 glider is a cycle running over four orbits. 
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Fig. 3. Glider over translation group T 2 is cycle in four group orbits. 



Our program constructs the maximum possible automorphism group for any 
lattice. For an N x N square toric lattice this group is the above mentioned 
Gnxn (we assume N 4, see formula (JTJ) and subsequent discussion). 

2 The program numbers orbits in the order of decreasing of their sizes and at equal 
sizes the lexicographic order of lexicograhically minimal orbit representatives is used. 

3 In Figs. [3] and [4] the configurations belonging to the same orbit have identical colors. 
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Now the glider is reproduced after two steps in time. As one can see from 
Fig. 01 ip^ is obtained from ipi and ^4 from tp2 by combinations of translations, 
90° rotations and reflections. Thus, the glider in torus (and in the discrete plane 
obtained from the torus as n — > 00) is a cycle located in two orbits of maximal 
automorphism group. 
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Fig. 4. Glider over maximal symmetry group T 2 X D4 is cycle in two group 
orbits. 



Note also that similar behavior is rather typical for continuous systems too. 
Many equations of mathematical physics have solutions in the form of running 
wave ip (x — vt) (= <p (<? t _1 :c) for Galilei group). One can see also an analogy 
between 11 spaceships" of cellular automata and solitons of KdV type equations. 
The solitons — like shape preserving moving structures in cellular automata — 
are often arise for rather arbitrary initial data. 

Cellular Automata with Symmetric Local Rules. As a specific class of 
discrete dynamical systems, we consider 'one-time-step' cellular automata on k- 
valent lattices with local rules symmetric with respect to all permutations of k 
outer vertices of the neighborhood. This symmetry property is an immediate dis- 
crete analog of general local diffcomorphism invariance of fundamental physical 
theories based on continuous space. The diffcomorphism group Diff(M) of the 
manifold M is very special subgroup of the infinite symmetric group Sym(M) 
of the set M. 

As we demonstrated in 0], in the binary case, i.e., if the number of vertex 
values q — 2, the automata with symmetric local rules are completely equivalent 
to generalized Conway's "Game of Life" automata [6] and, hence, their rules can 
be formulated in terms of "Birth"/ "Survival" lists. 

Adopting the convention that the outer points and the root point of the 
neighborhood are denoted x%, . . . , x% and Xk+i, respectively, we can write a local 
rule determining one-time-step evolution of the root in the form 

x' k+1 = f(x 1 ,...,x k ,x k+1 ). (4) 

The total number of rules (0| symmetric with respect to permutations of points 

/fc + q-l\ 

xx,-.. ,xt is equal to q\ 9- 1 > q . For the case of our interest (k = 3, q = 2) this 
number is 256. 
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It should be noted that the rules obtained from each other by permutation of 
q elements in the set Q are equivalent since such permutation means nothing but 
renaming of values. Thus, we can reduce the number of rules to consider. The 
reduced number can be counted via Burnside 's lemma as a number of orbits of 
rules ([4]) under the action of the group S q . The concrete expression depends on 
the cyclic structure of elements of S q . For the case q — 2 this gives the following 
number of non-equivalent rules 

N rules = 2 2fc+1 + 2 fe . 

Thus, studying 3-valent binary case, we have to consider 136 different rules. 

Example of Phase Portrait. Cellular Automaton 86. As an example 
consider the rule 86 on hexahedron. The number 86 is the "little endian" repre- 
sentation of the bit string 01101010 taken from the last column of the rule table 
with S3-symmetric combinations of values for xi,X2,X3 



X\ X 2 X 3 X4 


X4 








1 


1 


10 


1 


10 1 





110 


1 


110 1 





1110 


1 


1111 






The rule can also be represented in the "Birth"/ "Survival" notation as B123/S0, 
or as polynomial over the Galois field F 2 (see [5]) 

x' 4 = X4, + er 3 + er 2 + <J\ , 

where a\ — x\ + x 2 + X3, cr 2 = x\x 2 + xix 3 + X2X3, 03 = X1X2X3 are symmet- 
ric functions. In Fig. [5] the group orbits are represented by circles. The ordinal 
numbers of orbits are placed within these circles. The numbers over orbits and 
within cycles are sizes of the orbits (recall that all orbits included in one cycle 
have the same size). The rational number p indicates the weight of the corre- 
sponding element of phase portrait. In other words, p is a probability to be in 
an isolated cycle or to be caught by an attractor at random choice of state: p 
= (size of basin) / (total number of states). Here size of basin is sum of sizes of 
orbits involved in the struture. 

Note that most of cycles in Fig. [5] (36 of 45 or 80%) are "spaceships" . Other 
computed examples also confirm that soliton-like moving structures are typical 
for cellular automata. 

Of course, in the case of large lattices it is impractical to output full phase 
portraits (the program easily computes tasks with up to hundreds thousands of 
different structures). But it is not difficult to extract structures of interest, e.g., 
"spaceships" or "gardens of Eden" . 
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Fig. 5. Rule 86. Equivalence classes of trajectories on hexahedron. 36 of 45 cycles 
are 11 spaceships" . 

Search for Reversibility. The program is able to select automata with prop- 
erties specified at input. One of such important properties is reversibility. 

In this connection we would like to mention recent works of G. 't Hooft. One of 
the difficulties of Quantum Gravity is a conflict between irreversibility of Gravity 

— information loss (dissipation) at the black hole horizon — with reversibility 
and unitarity of the standard Quantum Mechanics. In several papers of recent 
years (see, e.g., |7|8j ) 't Hooft developed the approach aiming to reconcile both 
theories. The approach is based on the following assumptions 

— physical systems have discrete degrees of freedom at tiny (Planck) distance 
scales; 

— the states of these degrees of freedom form primordial basis of Hilbert space 
(with nonunitary evolution); 

— primordial states form equivalence classes: two states are equivalent if they 
evolve into the same state after some lapse of time; 

— the equivalence classes by construction form basis of Hilbert space with uni- 
tary evolution described by time-reversible Schrodinger equation. 

In our terminology this corresponds to transition to limit cycles: in a finite time 
of evolution the limit cycle becomes physically indistinguishable from reversible 
isolated cycle — the system "forgets" its pre-cycle history. Fig. illustrates 
construction of unitary Hilbert space from primordial. 

This irreversibility hardly can be found experimentally (assuming, of course, 
that considered models can be applied to physical reality). The system should 
probably spend time of order the Planck one (« 10~ 44 sec) out of a cycle and 
potentially infinite time on the cycle. Nowadays, the shortest experimentally 
fixed time is about 10 -18 sec or 10 26 Planck units only. 

Applying our program to all 136 symmetric 3-valent automata we have the 
following. There are two rules trivially reversible on all lattices 
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Primordial basis 
ei,e 2 ,e3,e4,e 5 ,e 6 ,e 7 



Equivalence classes 
Ei = {e 1 ,e 5 ,e 6 ,e 7 } 
E 2 = {e 2 } 




E 3 = {e 3 ,e 4 } 
form unitary basis 




Fig. 6. Transition from primordial to unitary basis. 



- 85 - B0123/S ~ 4 = x 4 + 1, 

- 170 ~ B/S0123 ~ x' 4 = x 4 . 

Besides these uninteresting rules there are 6 reversible rules on tetrahedron 



- 43 - 


- B0/S012 r 


° 


= x 4 (a 2 4 


-^)- 


f cr 3 + (72 - 


+- o-i + 1, 


- 51 - 


- B02/S02 r 


^4 


= o-i + l, 








- 77 - 


- B013/S1 r 


° 3^4 


= x 4 (<7 2 4 


■ o~i -1 


- 1) + cr 3 + 


O 2 + 1, 


- 178 r 


- B2/S023 - 


^4 


= ^4 (^2 - 


-CTl H 


- 1) + cr 3 4 


- er 2 , 


- 204 r 


- B13/S13 ' 


^4 










- 212 r 


- B123/S3 - 


X4 


= x 4 (a 2 - 


'rax) 


+ (T3 + (72 





Note that all these reversible rules are symmetric with respect to permutation 
of values Q — {0, 1}. Two of the above rules, namely 51 and 204, are reversible 
on hexahedron too. There are no nontrivial reversible rules on all other lattices 
from Fig. Q] Thus we may suppose that 't Hooft's picture is typical for discrete 
dynamical systems. 

5 Statistical Lattice Models and Mesoscopic Systems 

Statistical Mechanics. The state of deterministic dynamical system at any 
point of time is determined uniquely by previous states of the system. A Markov 
chain — for which transition from any state to any other state is possible with 
some probability — is a typical example of non- deterministic dynamical system. 
In this section we apply symmetry approach to the lattice models in statistical 
mechanics. These models can be regarded as special instances of Markov chains. 
Stationary distributions of these Markov chains are studied by the methods of 
statistical mechanics. 

The main tool of conventional statistical mechanics is the Gibbs canonical 
ensemble - imaginary collection of identical systems placed in a huge thermostat 
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with temperature T. The statistical properties of canonical ensemble are encoded 
in the canonical partition function 

Z=J2 c- E ° /kBT ■ (5) 

Here S is the set of microstates, E a is energy of microstate a, ks is Boltzmann's 
constant. The canonical ensemble is essentially asymptotic concept: its formula- 
tion is based on approximation called "thermodynamic limit" . For this reason, 
the canonical ensemble approach is applicable only to large (strictly speaking, 
infinite) homogeneous systems. 

Mesoscopy. Nowadays much attention is paid to study systems which are too 
large for a detailed microscopic description but too small for essential features 
of their behavior to be expressed in terms of classical thermodynamics. This 
discipline, often called mesoscopy, covers wide range of applications from nuclei, 
atomic clusters, nanotechnological structures to multi-star systems |9|10lllj . To 
study mesoscopic systems one should use more fundamental microcanonical en- 
semble instead of canonical one. A microcanonical ensemble is a collection of 
identical isolated systems at fixed energy. Its definition does not include any 
approximating assumptions. In fact, the only key assumption of a microcanon- 
ical ensemble is that all its microstates are equally probable. This leads to the 
entropy formula 

S E = k B In fl E , (6) 
or, equivalently, to the microcanonical partition function 

Q E = c S E /k B _ (y) 

Here Qe is the number of microstates at fixed energy E. In what follows we 
will omit Boltzmann's constant assuming k B = 1. Note that in the thermody- 
namic limit the microcanonical and canonical descriptions are equivalent and 
the link between them is provided by the Laplace transform. On the other hand, 
mesoscopic systems demonstrate observable experimentally and in computation 
peculiarities of behavior like heat flows from cold to hot, negative specific heat 
or "convex intruders" in the entropy versus energy diagram, etc. These anoma- 
lous - from the point of view canonical thermostatistics - features have natural 
explanation within microcanonical statistical mechanics |llj . 

Lattice Models. In this section we apply symmetry analysis to study meso- 
scopic lattice models. Our approach is based on exact enumeration of group 
orbits of microstates. Since statistical studies are based essentially on different 
simplifying assumptions, it is important to control these assumptions by exact 
computation, wherever possible. Moreover, we might hope to reveal with the help 
of exact computation subtle details of behavior of system under consideration. 
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As an example, let us consider the Ising model. The model consists of spins 
placed on a lattice. The set of vertex values is Q = { — 1,1} and the interaction 
Hamiltonian is given by 

H = -J^8i8j-B^8i , (8) 

(ij) ' 

where Sj, Sj ; 6 Q; J is a coupling constant (J > and J < correspond to 
ferromagnetic and antiferromagnetic cases, respectively); the first sum runs over 
all edges (i,j) of the lattice; B is an external "magnetic" field. The second sum 
M = Si is called the magnetization. To avoid unnecessary technical details we 

i 

will consider only the case J > (assuming J = 1) and B = in what follows. 

Since Hamiltonian and magnetization are constants on the group orbits, we 
can count numbers of microstates corresponding to particular values of these 
functions - and hence compute all needed statistical characteristics - simply by 
summation of sizes of appropriate orbits. 

Fig. [7] shows microcanonical partition function for the Ising model on do- 
decahedron. Here total number of microstates fl — 1048576, number of lattice 
vertices = 20, energy E is value of Hamiltonian. 

0.15 | , , , , , , , , , 1 

0.125 - 

■ 

0.1 - 

_ 0.075 - 

0.05 - • 

0.025 - 

I . 1 ■ 1 . 1 . 1 .-! . 1 

-1.5 -1 -0.5 0.5 1 

e 

Fig. 7. Microcanonical density of states p(e) = QeI ^ versus energy per vertex 
e = EJ \ V(r)\ for the Ising model on dodecahedron. 



Of course, other characteristics of the system can be computed easily via 
counting sizes of group orbits. For example, the magnetization is shown in Fig. 

m 
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Fig. 8. Specific magnetization 771(e) = M(E)/ \V(r)\ vs. energy per vertex e for 
the Ising model on dodecahedron. 

Phase Transitions. Needs of nanotechnological science and nuclear physics 
attract special attention to phase transitions in finite systems. Unfortunately 
classical thermodynamics and the rigorous theory of critical phenomena in ho- 
mogeneous infinite systems fails at the mesoscopic level. Several approaches have 
been proposed to identify phase transitions in mesoscopic systems. Most accepted 
of them is search of "convex intruders" [12] in the entropy versus energy diagram. 
In the standard thermodynamics there is a relation 



where Cy is the specific heat at constant volume. It follows from ([9]) that 
d 2 S/dE 2 \ < and hence the entropy versus energy diagram must be concave. 
Nevertheless, in mesoscopic systems there might be intervals of energy where 
d 2 S/dE 2 \ v > 0. These intervals correspond to first-order phase transitions and 
are called "convex intruders". From the point of view of standard thermodynam- 
ics one can say about phenomenon of negative heat capacity, of course, if one 
accepts that it makes sense to define the variables T and Cy as temperature and 
the specific heat at these circumstances. In [13] it was demonstrated via compu- 
tation with exactly solvable lattice models that the convex intruders flatten and 
disappear in the models with local interactions as the lattice size grows, while 
in the case of long-range interaction these peculiarities survive even in the limit 
of an infinite system (both finite and long-range interacting infinite systems are 
typical cases of systems called nonextensive in statistical mechanics). 



d 2 s 



1 1 



(9) 



BE 2 



v 



14 



A convex intruder can be found easily by computer for the discrete systems 
we discuss here. Let us consider three adjacent values of energy Ei-i,Ei,Ei + i 
and corresponding numbers of microstates J7e 4 _u flEn ^E i+1 - In our discrete 
case the ratio (-E/j+i — Ei) / (Ei — Ei-\) is always rational number p/q and we 
can write the convexity condition for entropy in terms of numbers of microstates 
as easily computed inequality 

rt£'<^_A +1 • (io) 

As a rule Ei + i — Ei = Ei — Ei-\ and inequality (fTU)) takes the form 

o% < f2 El _ 1 n Ei+1 . 

This form means that within convex intruder the number of states with the 
energy Ei is less than geometric mean of numbers of states at the neighboring 
energy levels. 

Fig. [5] shows the entropy vs. energy diagram for the Ising model on do- 
decahedron. The diagram has apparent convex intruder in the energy interval 
[—24, —18]. Exact computation reveals also a subtle convex intruder in the in- 
terval [—16, —12]. (In terms of specific energy, as in Fig. [51 these intervals are 
[—1.2, —0.9] and [—0.8, —0.6], respectively.) It is well known that one-dimensional 
Ising model has no phase transitions. To illustrate the difference between the di- 
agrams for the cases with and without phase transitions, we place also in Fig. [9] 
the diagram for Ising model on the ID circle lattice with 24 vertices. 
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Fig. 9. Specific microcanonical entropy s(e) = ln(f?E) / |V^(.T)| vs. energy per 
vertex e for the Ising model on dodecahedron (left) and on circle of length 24 
(right). Left diagram contains distinct convex intruder in the interval —1.2 < 
e < —0.9 and subtle one in the interval —0.8 < e < —0.6. Right diagram is fully 
concave: one-dimensional Ising model has no phase transitions. 

In Fig. [10] we show the entropy-energy diagrams for lattices of different va- 
lences, namely, for 3-, 4- and 6-valent tori. These lattices are marked in Fig. 
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Q]as "Graphene 6x4", "Square 5x5" and "Triangular 4x6", respectively. The 
diagram for 3-valent torus is symmetric with respect to change sign of energy 
and contains two pairs of adjacent convex intruders. One pair lies in the e- 
interval [—1.25,-0.75] and another pair lies symmetrically in [0.75,1.25]. The 
4-valent torus diagram contains two intersecting convex intruders in the inter- 
vals [—1.68,-1.36] and [—1.36,-1.04]. The 6-valent torus diagram contains a 
whole cascade of 5 intersecting or adjacent intruders. Their common interval is 
[-2.5,-0.5]. 
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Fig. 10. Specific microcanonical entropy for the Ising model on 3-valent (dot 
line, 24 vertices), 4-valent (dash line, 25 vertices) and 6-valent (solid line, 24 
vertices) tori. 



6 Summary 

— AC program for symmetry analysis of finite discrete dynamical systems has 
been created. 

— We pointed out that trajectories of any deterministic dynamical system go 
always in the direction of nondecreasing sizes of group orbits. Cyclic trajec- 
tories run within orbits of the same size. 

— After finite time evolution operators of dynamical system can be reduced to 
group actions. This lead to formation of moving soliton-like structures - 
"spaceships" in the case of cellular automata. Computer experiments show 
that "spaceships" are typical for cellular automata. 
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— Computational results for cellular automata with symmetric local rules allow 
to suppose that reversibility is rare property for discrete dynamical systems, 
and reversible systems are trivial. 

— We demonstrated capability of exact computing based on symmetries in 
search of phase transitions for mesoscopic models in statistical mechanics. 
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